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Abstract 

Electronic properties of graphene oxides enriched by the strong chemical bondings 
are investigated using first-principle calculations. They are very sensitive to the 
changes in the number of graphene layer, stacking configuration, and distribution of 
oxygen. The feature-rich electronic structures exhibit the destruction or distortion 
of Dirac cone, opening of band gap, anisotropic energy dispersions, O- and (C,0)- 
dominated energy dispersions, and extra critical points. All the few-layer graphene 
oxides are semi-metals except for the semiconducting monolayer ones. For the former, 
the distorted Dirac-cone structures and the O-dominated energy bands near the Fermi 
level are revealed simultaneously. The orbital-projected density of states (DOS) have 
many special structures mainly coming from a composite energy band, the parabolic 
and partially flat ones. The DOS and spatial charge distributions clearly indicate the 
critical bondings in O-O, C-0 and C-C bonds, being responsible for the diversified 
properties. 
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1. Introduction 


Graphene has been a mainstream material in science and engineering for so many years, 
mainly owing to its remarkable physical, chemical and material properties. Numerous re¬ 
searchers have tried to diversify the essential properties of graphene by using doping M, 
stacking configuration [HE], layer number [BJC], electric or/and magnetic fields [SHIT], and 
mechanical strain [T21H3] . Monolayer graphene has a low-energy isotropic Dirac-cone struc¬ 
ture, while it is a zero-gap semiconductor because of the vanishing density of states (DOS) 
at the Fermi level (Ep= 0). On the other hand, all the few-layer graphenes are semi-metals 
with small overlaps in valence and conduction bands, in which the low-lying energy bands 
are significantly affected by the layer number and stacking configurations. Recently, oxygen 
absorbed graphene as so-called graphene oxide (GO) has been the focus of attention due to 
its opening of band gap and other interesting properties [HHl6j. It has promised potential 
applications in wide variety of areas, including memristor devices mi. sensors P2H32], and 
supercapacitors |20F[22| . The paper seeks to contribute to critical orbital hybridizations 
in O-O, C-0 and C-C bonds by investigating how oxygen absorption, layer number and 
stacking configuration can enrich the electronic properties of graphene. 

The well known method to produce GO is the Hummers method in which a mixture of 
sodium nitrate, sulphuric acid and potassium permanganate are used in the synthesis pro¬ 
cess [23] , Recently, other methods have been developed by adding phosphoric acid combined 
with the sulphuric acid without sodium nitrate, and increasing the amount of potassium 
permanganate [24]. The oxygen concentration of 70% can be reached in this modified 
method [24]. The distribution and concentration of oxygen can be examined using nuclear 
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magnetic resonance (NMR) [25], and X-ray photoelectron spectroscopy (XPS) [55], respec¬ 
tively. On the other hand, theoretical studies have shown the dependence of band gap on 
the concentration and distribution of oxygen atoms in monolayer graphene oxide [271 [28] . 
The oxygen absoption has converted a part of sp 2 carbon bonds in pristine graphene to 
sp 3 bonds in GO [35]. It has been reported that the strong hybridization between the 2 p x>y 
orbitals of O atoms and the 2 p z orbitals of C atoms is responsible for the opening of band 
gap [57] . However, this work shows that the above-mentioned orbital hybridization is not 
exact (discussed later in Figs. 4 and 5 ). The important orbital hybridizations in O-O, C-0 
and C-C bonds are not fully analyzed, and the oxygen absorbed few-layer graphenes sys¬ 
tem has not been investigated. Therefore, the critical roles of such bondings in determining 
the essential properties of oxygen absorbed few-layer graphenes are worthy of a systematic 
study. 

In this paper, the strong effects of oxygen atoms on the electronic properties of graphene, 
with different distributions, layers, and stackings, are investigated using first-principles cal¬ 
culations. The calculated results demonstrate that the bond lengths, energy dispersions, 
band gap, charge distributions, and density of states (DOS) are sensitive to the changes in 
oxygen distributions and geometric structures. Specifically, the energy bands exhibit a lot 
of features: the destruction or distortion of Dirac cone, opening of band gap, partially flat 
bands, energy dispersions related to 0-0 and C-0 bonds, and extra critical points. Such 
features are reflected in the orbital-projected DOS, including the absence and presence 
of 7r band peaks, the O-dominated special structures near E F , and the (C,0)-dominated 
prominent peaks at middle energy. Furthermore, the diverse orbital hybridizations in C-0, 
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0-0, and C-C bonds will be explored through the DOS and spatial charge distributions 
to explain the dramatic changes in electronic properties. The above-mentioned results can 
be verified by experimental measurements such as ARPES [29], and scanning tunneling 
spectroscopy (STS) [7]. 


2. Methods 

The first-principles calculations on GO are performed based on density functional the¬ 
ory (DFT) using the Vienna Ab initio Simulation Package (VASP) [3011.31] . The electron- 
ion interactions are evaluated by the projector augmented wave method [32], whereas the 
electron-electron interactions are taken into account by the exchange-correlation function 
under the generalized gradient approximation of Perdew-Burke-Ernzerhof [33]. A vacuum 
layer with a thickness of 12 A is added in a direction perpendicular to the GO plane to avoid 
interactions between adjacent unit cells. The wave functions are expanded using a plane- 
wave basis set with a maximum kinetic energy of 500 eV. All atomic coordinates are relaxed 
until the Hellmann-Feynman force on each atom is less than 0.01 eV/A. Furthermore, the 
van der Waals (vdW) force is employed in the calculations using the semiemprical DFT- 
D2 correction of Grimme to correctly describe the atomic interactions between graphene 
layers [34] , The k-point mesh is set as 200x200x1 in geometry optimization, 30x30x1 in 
band structure, and 250x250x1 in the DOS for the lxl unit cell with zigzag structure. 
An equivalent k-point mesh is set for other cells depending on their size. 


3. Results and discussion 
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The atomic structures of oxygen absorbed monolayer graphene with top view and side 
view are shown in Fig. 1. Oxygen atoms are absorbed on the top graphene layer with 
O/C ratio of 1/2 (50%) for zigzag (Fig. 1(a)), armchair (Fig. 1(b)), and chiral (Fig. 1(c)) 
unit cells, in which carbon atoms are arranged along these edge structures, respectively. 
In comparison with the two latters, the first one has a lower total ground state energy 
meaning more stable. This is in agreement with previous studies which have shown that 
the higher the symmetric of the geometric structure, the more stable the system is [35] . 
However, the difference in total ground state energy among these configurations is only 
about 7-8 %, indicating that all the distributions are considerably stable. From the top 
view, oxygen atoms are at the top between two carbon atoms as so-called the bridge-site. 
We have checked the most stable sites when oxygen is absorbed on graphene by comparing 
their total ground state energies. Among the bridge-, hollow- and top-sites, the bridge-site 
is the most stable one and the hollow-site is the less stable one, which agrees with previous 
studies 

The main features of geometric structures, including the C-C bond lengths, C-0 bond 
lengths, and interlayer distances are dominated by the number of graphene layer (n). As 
oxygen atoms are absorbed on graphene, the C-C bond lengths are expanded compared 
to 1.424 A in pristine graphene, as shown in Table 1. The equivalent C-C bond lengths 
are altered after oxidation, in which the nearest C-C bond corresponding to oxygen atom 
on the bridge side is shorter than the second nearest one, as a result of the strong C-0 
bond (Fig. 4(b)). The longer C-C bond length in GO means the 2 p x , y orbitals have the 
hybridization with atomic orbitals of O. As the layer number grows n>2, the C-C bond 


5 


lengths show a significant decrease while the C-0 bond length only has a slight rise. They 
remain almost unchanged and close to that of pristine graphene for n>8. Our calculation 
also reveals that the interlayer distances are affected by oxygen absorption. The optimized 
values for AA and AB stackings are equal to 3.27 A instead of 3.52 A and 3.26 A in pristine 
bilayer systems [3] . Similarly, for AAA and ABA stackings, the interlayer distance between 
the middle and top layers is about 3.27 A, whereas that of the bottom and middle layers 
are nearly identical to the AA- and AB-bilayer ones, respectively. Further examinations, 
it is found that the n-layer graphene with oxygen absorbed on the top can be qualitatively 
regarded as the (n—l)-layer of graphene with monolayer graphene oxide. The chemical 
bonding of atomic orbitals will be investigated in more detail later (Fig. 4). 

The two-dimensional (2D) energy bands along high symmetric points are determined 
by the C-O, 0-0 and C-C bonds, oxygen distribution, layer number and stacking config¬ 
uration (Fig. 2). The contributions of the O atoms and the C atoms passivated with the 
former are represented by the blue and red circles, respectively, in which the dominance 
is proportional to the circle’s radius. Different from monolayer graphene, the isotropic 
Dirac-cone structure near the K point is destroyed in GO, mainly owing to the serious 
hybridization between the atomic orbitals of C and O atoms (discussed later in charge 
density and DOS). Such strong C-0 bonds lead to the termination of the complete 7r bonds 
between parallel 2 p z orbitals of C atoms accounting for the Dirac-cone structure near Ep in 
pristine graphene. Instead, there is a wide direct gap of 3.5 eV, 2.8 eV, and 3.9 eV for the 
zigzag, armchair, and chiral unit cells, as shown in Figs. 2(a), 2(c) and 2(e), respectively. 
The highest occupied state is related to an O-dominated energy band. The 0-0 bond 
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has the weaker a bond which arises from the (2p y —2p y ), (2p xy —2p Xjy ) and {2p x ^—2p xy ) 
hybridizations for zigzag, armchair and chiral distributions, respectively (Figs. 4 and 5). 
As a result, the O-dominated low-lying bands might have the partially flat or parabolic 
dispersions for the different distributions. On the other hand, the above-mentioned C-0 
bonds induce the (C,0)-related energy bands at -2 eV< E v < -4 eV because of the strong 
orbital hybridizations between passivated C and O atoms (blue and red circles, revealed si¬ 
multaneously). As to the planar C-C bond, it possesses the strongest a bond formed by the 
(2 p Xl 2p yi 2s) orbitals and thus creates the deeper a bands with E v < —3.5 eV (red circles). 

For bilayer GO, the band gap is replaced by a distorted Dirac-cone structure that comes 
from the C atoms on another layer without O passivation. The Dirac point is revealed at 
the K point for zigzag (Fig. 2(b)) and chiral unit cells (Fig. 2(f)) and at the T point 
for armchair one (Fig. 2(d)). The anisotropic linear bands will change into the parabolic 
bands with the saddle point along K —)■ M or T —> M. For example, the saddle points of 
the 7r and n* bands, respectively, correspond to E v « —2.5 eV and E c « 1.3 eV for the 
zigzag distribution in Fig. 2(b). Different from AA bilayer graphene in which two isotropic 
Dirac-cone structures with the linear bands intersecting at Ep [5], there is one anisotropic 
Dirac-cone structure with the Fermi level located at the conduction band. There still exist 
the O-dominated energy dispersions near Ep and the (C,0)-related bands at middle en¬ 
ergy, as revealed in monolayer GO. The band-edge states of the former are above the Fermi 
level, indicating its free-hole density equal to the free-electron one in the distorted Dirac 
cone. As the number of graphene layer increases, there are more distorted Dirac cones and 
saddle points. The band structure differences between AA and AB stackings (Fig. 2(g)) 
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are small and only lie in the crossing and anti-crossing bands (green ellipses). Specially, 
the energy dispersions of ABA and ABC stackings (Fig. 2(k)) are almost identical, since 
these two systems correspond to the same structure of AB bilayer graphene with monolayer 
GO on the top, as proposed ealier in the geometric structures. However, these two systems 
are different from the AAA stacking in the distorted Dirac-cone structures as well as the 
crossing and anti-crossing bands (Figs. 2(h) and 2(i)). In addition, the Fermi velocities of 
the distorted Dirac cones are reduced about 7-10% compared with that of pristine graphene 
(» 10 6 m/s) j37]. 

The critical low-energy electronic properties can be further understood by the three- 
dimensional (3D) bands shown in Fig. 3. The zigzag AA stacking, respectively, exhibits 
the partially flat and parabolic dispersions on the xz and yz planes (blue and green re¬ 
gions) with the crossing of the Fermi level. Therefore, the O-dominated energy band can 
be regarded as a ID parabolic band. However, this band has the parabolic (the partially 
flat dispersions) on the planar projections for the armchair (chiral) distribution, as shown 
in Fig. 2(c) (Fig. 2(e)). As to the distorted Dirac-cone structures, their main features 
include the change in energy dispersion, the separation of Dirac points, and the crossing 
with the O-dominated energy band. For example, the armchair AA stacking has an energy 
spacing of 50 meV between two separate Dirac points (near the F point in Fig. 3(b)) and 
one pair of O-dominated parabolic bands at E v ~ —0.75 eV crossing the valence Dirac 
cone (Fig. 2(d)). The low-lying energy bands of zigzag AAA and ABA stackings, as shown 
in Figs. 3(c) and 3(d), clearly illustrate the significant differences in energy dispersions 
and spacings. The former are two pairs of linear bands with ~ 10 meV energy spacing, 


whereas the later consist of one pair of parabolic bands and one pair of Mexican-hat bands 
with ~ 0.1 eV energy spacing. 

Recently, ARPES has emerged as a most useful experimental technique to study the 
electronic band structures. The experimental measurements on graphene-related systems 
have been used to investigate the effects due to doping [38], layer number stacking 

configuration [29], and the electric held [29] . For example, they have verified the Dirac-cone 
structure of graphene grown on SiC [10] , and observed the opening of band gap for graphene 
on fr(lll) through oxidation [41] . As expected, the feature-rich energy bands of oxygen 
absorbed few-layer graphenes, including the absence and presence of the distorted Dirac- 
cone structures, the band gap, the O-dominated bands near E F , and the (C,0)-dominated 
bands at middle energy can be examined by ARPES. The comparisons between theoreti¬ 
cal predictions and experimental measurements can comprehend how oxygen distribution, 
layer number and stacking configuration affect the electronic properties of oxygen absorbed 
few-layer graphenes. 

The charge density (p) and the charge density difference (Ap) can provide very useful 
information on the spatial charge redistributions and thus on the dramatic changes of en¬ 
ergy bands [42] - The former reveals the chemical bondings as well as the charge transfer. 
In pristine A A stacking (Fig. 4(a)), a strong covalent a bond with high charge density 
due to the ( 2p x ,2p y ,2s ) orbitals exists between two C atoms (the pink rectangle; details in 
Fig. 4(d)). Furthermore, there is a weak tt bond arising from the parallel 2 p z orbitals at 
the boundary (the purple rectangle). When O atoms are absorbed on zigzag AA stacking 
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surface (Fig. 4(b)), charges transferred from C to O atoms are ~ le using the Bader 
charge analysis (the dashed black rectangle). The strong C-0 bond, respectively, induces 
the termination of 7 r bond on the upper boundary and the charge redistribution of n bond 
on another one. The former accounts for the absence of Dirac cone and the opening of 
energy gap. The latter is further reflected in the 7r bond of the bottom layer by the vdW 
interactions, so that the Dirac cone presents the distorted structure (Fig. 2(b)). The 
above-mentioned characteristics about the o and 7r bonds as well as the strong C-0 bond 
are also found in armchair and chiral distributions (not shown). 

In order to comprehend the orbital hybridizations in O-O, C-0 and C-C bonds, the 
variation of charge density of zigzag AA bilayer is illustrated in Figs. 4(c) and 4(d). A p 
is created by substracting the charge density of isolated C and O atoms from that of GO 
system. The top view of charge density difference of zigzag AA stacking is shown in the 
inset of Fig. 4(c), where the red and purple dashed lines represent the slices in yz and xz 
planes, respectively. The p y orbitals of O with high charge density enclosed in the dashed 
black rectangle (Fig. 4(c)) are lengthened along y compared to isolated O atom’s, clearly 
indicating the 2p y —2p y hybridization. The 0-0 bond belongs to a weak cr bond because 
of the larger distance, being the main reason for the low-lying O-related bands near Ep 
(Fig. 2(b)). It has also been observed that the 0-0 bond comes from the hybridization of 
^Px,y~^Px,y f° r the armchair and chiral distributions (not shown: Figs. 5(c) and 5(e)). In 
sharp contrast to the 2 p y orbitals, the 2 p z and 2 p x orbitals of O have strong hybridizations 
with those of C, as seen from the green region enclosed by the dashed purple rectangle 
(Fig. 4(d)). This region lies between O and passivated C atoms and bents toward the 
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latter, while the p y orbitals do not reveal any hybridization with other orbitals. In fact, it 
is hard to completely distinguish which orbitals have hybridizations to each other by using 
charge density analysis, but they also show the significant evidence in the orbital-projected 
DOS (Fig. 5). The C-0 bond is much stronger than the 0-0 bond; therefore, the (C,0)- 
related bands appear at the range of -2 eV< E v <-4 eV. On the other hand, between 
two non-passivated C atoms of GO, A p shows a strong a bond indicated by the enclosed 
pink square in Figs. 4(d). Such bond becomes a bit weaker after the C-0 bond is formed 
(the dashed pink square in Figs. 4(d)). This demonstrates that not only the 2 p z orbitals 
of passivated C atoms hybridize with orbitals of O atoms but also the 2 p x y orbitals play 
an important role. However, the planar C-C bond remains the strongest one and has the 
deeper a energy bands. 

The main characteristics of band structures are directly reflected in DOS. The orbital- 
projected DOS, as shown in Fig. 5, is useful in understanding the orbital contributions as 
well as the orbital hybridizations in chemical bonds. The feature-rich DOS presents three 
kinds of special structures, including the asymmetric peak (circle), the shoulder structures, 
and the symmetric peaks. They, respectively, correspond to a composite band with par¬ 
tially flat and parabolic dispersions, the minimum/maximum band-edge states of parabolic 
bands, and the saddle points of parabolic bands or the partially flat bands [12]. The low- 
energy DOS is thoroughly altered after oxygen absorption. For pristine graphene, the peaks 
caused by n and tt* bands due to 2p z -2p z bondings between C atoms will dominate within 
the range of \E C,V \ < 2 eV [37]. However, the tt- and 7r*-peaks are absent as a result of the 
strong C-0 bond. Instead, there are several O-dominated prominent structures in a wide 
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range of E ~-2.5 eV to Ep, mainly coming from the 2p y —2p y bondings in zigzag distri¬ 
bution (Fig. 5(a)) and the 2p Xiy —2p XiV bondings in armchair and chiral ones (Figs. 5(c) 
and 5(e)). It is noticed that the contributions of 2p x and 2 p y orbitals in armchair distribu¬ 
tion are equivalent because of the symmetric distribution of oxygen. With the addition of 
graphene layer, the it- and 7r*-peaks related to the non-passivated C atoms appear around 
-2.5 eV and 1.2 eV, respectively (arrows in Figs. 5(b), 5(d); 5(f)-5(h)). Apparently, the 
number of these peaks are only affected by graphene layer number (Figs. 5(b), 5(g) and 
5(h)). 

With the increasing energy, DOS grows quickly and exhibit prominent symmetric peaks 
due to the C-0 bond at middle energy (-2 eV< E v <-4 eV). Among these peaks, the most 
pronounced one is mainly contributed by the (2p z ,2p x ), (2p x ,2p y ,2p z ) and {2p xi 2p yi 2p z ) 
orbitals of C and O atoms for the zigzag, armchair and chiral distributions, respectively 
(e.g., triangles in Figs. 5(a), 5(c); 5(e)). The contributions of distinct orbitals revealed at 
the same energy apparently represent the significant hybridizations between them, being 
consistent with those in the spatial charge distributions (Fig. 4(d)). This strong covalent 
C-0 bond is the main reason for the absence of tt and tt* peaks and the opening of band 
gap. Also, the magnitude of energy gap depends on the 0-0 bond strength, as indicated 
from distinct gaps among three types of distributions. 

Up to now, there are several theoretical studied on GO, especially in electronic proper¬ 
ties [E3IZ3I2B1II3]- They are focused on how the band gap is modulated with the variation 
of oxygen concentrations I2ZII2S]. The orbital hybridizations that cause the opening of band 
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gap and critical electronic properties have not been fully explored. In this work, the com¬ 
prehensive orbital hybridizations in C-O, O-O, and C-C bonds are analyzed by the spatial 
charge distributions and orbital-projected DOS (Figs. 4 and 5). The previous study [27] 
proposed 2p z —2p xy orbital hybridizations of C and O atoms to explain the opening of 
band gap, in which they are not totally exact. Our results indicate orbital hybridizations 
of 2p X}Z —2p xz or 2p x y z —2p Xt y Z between passivated C and O atoms, obviously including 
the 2p z -2p z hybridizations. The effect of oxygen distributions on the magnitude of band 
gap has been pointed out [U,;2T], but not the chemical picture for this, the 2p y —2p y or 
2Px,y—2Px,y hybridizations in 0-0 bond. Moreover, the terminated tt bond, the reformed 
7 r and a bonds are clearly illustrated using spatial charge distributions (Fig. 4). 

The STS measurement, in which the tunneling conductance (dl/dV) is approximately 
proportional to DOS and directly reflects its special structures, can provide an efficient 
way to confirm the distribution of oxygen. This method has been used for adatoms on 
graphene [44], carbon nanotubes [45], graphene nanoribbons [46], and few-layer graphene 
m For example, it has been used to observe the Fermi-level shift and additional peaks 
near Dirac point of graphene irradiated with Ar + ions j48j. The main features in electronic 
properties, including the energy gaps, the O-dominated prominent structures near Ep, the 
high 7r- and 7r*-peaks, and the strong (C,0)-related peaks at middle energy, can be further 
investigated with STS. The STS measurements on the low- and middle-energy peaks can 
identify the complex chemical bondings in oxygen absorbed few-layer graphenes. 


4. Conclusion 
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In summary, the geometric structures and electronic properties of oxygen absorbed few- 
layer graphenes are studied using first-principles calculations. They are shown to be dom¬ 
inated by the diverse chemical bondings between atomic orbitals. The C-C bond lengths 
are expanded due to oxygen absorption and gradually recover to that in pristine graphene 
with the increasing layer number. The nearest C-C bond length is shorter compared to the 
second nearest one, indicating the significant interactions between the passivated C atoms 
and O atoms. The C-C, C-0 and 0-0 bonds are responsible for the dramatic changes 
in electronic properties, including the destruction or distortion of Dirac cone, energy gap, 
anisotropic energy spectra, (C,0)- and O-dominated energy bands, and many extra critical 
points. The oxygen absorbed few-layer graphenes, with the O-related energy bands and 
the distorted Dirac cone near Ep, are semi-metals except for the semiconducting monolayer 
ones. The predicted rich-feature band structures, which depend on the oxygen distribution, 
staking configuration and layer number, could be examined by the ARPES measurements. 

The competition or cooperation among the critical chemical bondings in C-C, 0-0 
and C-0 bonds can enrich the essential properties. The complete 7r bondings due to the 
parallel 2p z orbitals in C-C bonds can be formed on the graphene planes without oxygen 
adsorption. They are affected by the oxygen passivation so that the distorted Dirac-cone 
structures exhibit the crossing with the O-dominated band, the energy spacings between 
the separated Dirac points, and the changes in energy dispersions. Such structures are 
sensitive to the number of layers and stacking configurations, e.g., the same bands between 
ABA and ABC stackings, but the distinct bands between ABA and AAA stackings. In 
addition, the (2p x ,2p y ,2s) orbitals of C atoms only create the a bands with the deeper 
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state energies lower than —3.5 eV. The 2p y —2p y or 2p X)y —2p Xty orbital hybridizations in 
the 0-0 bonds are determined by the oxygen distribution; furthermore, the sufficiently 
strong bondings will generate the O-dominated energy dispersions with a bandwidth of 
~2.5 eV near Ep. There exist serious orbital hybridizations of 2p xz —2p xz or 2p xiV ^ z —2p x V:Z 
between passivated C and O atoms, thus, leading to the deeper (C,0)-related energy bands 
at -2 eV< E v <-4 eV, and the absence of the isotropic Dirac cone, which has not been 
sufficiently explored in previous studies (T5ll2Tj . The main features of the orbital-dependent 
energy bands are directly reflected in a lot of peaks and shoulder structures in DOS. The 
experimental examinations of STS on DOS, as well as those of ARPES on energy bands, 
can provide the full information to understand the complex orbital bondings in graphene 
oxides. 
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FIGURE CAPTIONS 


Fig. 1. Geometric structures of 50% oxygen concentration systems with: (a) zigzag, 
(b) armchair, and (c) chiral unit cells. 

Fig. 2. Band structures of the oxygen-absorbed few-layer graphenes: (a) monolayer 
(zigzag), (b) AA bilayer (zigzag), (c) monolayer (armchair), (d) AA bilayer (armchair), (e) 
AB bilayer (zigzag), (f) AAA trilayer (zigzag), (g) ABA trilayer (zigzag), and (h) ABC 
trilayer (zigzag). Superscripts c and v correspond to the conduction and valence bands, 
respectively. Also shown in the inset of (a) is the first Brillouin zone. 

Fig. 3. Low-lying 3D band structures for: (a) zigzag AA bilayer near the M point, (b) 
armchair AA bilayer near the Gamma point, (c) zigzag AAA trilayer near the K point, and 
(d) zigzag ABA trilayer near the K point. 

Fig. 4. The charge density p of (a) pristine AA stacking, and (b) zizag AA bilayer. The 
charge density difference A p of zigzag AA bilayer in (c) yz-plane, and (d) xz-plane. 

Fig. 5. DOS of the oxygen-absorbed few-layer graphenes: (a) zigzag monolayer, (b) 
zigzag AA bilayer, (c) armchair AA bilayer, (d) chiral AA bilayer, (e) zigzag AB bilayer, 
and (f) zigzag AAA trilayer. 
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Figure 1: Geometric structures of 50% oxygen concentration systems with: (a) zigzag, (b) 
armchair, (c) chiral unit cells. 
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Table 1: The calculated C-0 bond lengths, and C-C bond lengths of the zigzag, armchair, 


and chiral unit cell with various layers. 


Graphene Oxides 

C-0 bond 

Nearest C-C bond 

Second nearest C-C 


length (A) 

length (A) 

bond length (A) 

Monolayer (zigzag) 

1.434 

1.503 

1.543 

Bilayer 

1.436 

1.459 

1.491 

Trilayer 

1.437 

1.441 

1.470 

Four-layer 

1.438 

1.432 

1.459 

Five-layer 

1.438 

1.431 

1.453 

Six-layer 

1.438 

1.430 

1.450 

Seven-layer 

1.438 

1.423 

1.449 

Eight-layer 

1.439 

1.418 

1.444 

Nice-layer 

1.439 

1.416 

1.444 

Monolayer (armchair) 

1.410 

1.543 

1.569 

Bilayer (armchair) 

1.413 

1.477 

1.497 

Monolayer (chiral) 

1.424 

1.519 

1.560 

Bilayer (chiral) 

1.429 

1.465 

1.496 
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Figure 2: Band structures of the oxygen-absorbed few-layer graphenes: (a) monolayer 
(zigzag), (b) AA bilayer (zigzag), (c) monolayer (armchair), (d) AA bilayer (armchair), (e) 
AB bilayer (zigzag), (f) AAA trilayer (zigzag), (g) ABA trilayer (zigzag), and (h) ABC 
trilayer (zigzag). Superscripts ”c” and ”v” correspond to the conduction and valence band, 
respectively. Also shown in the inset of (a) is the first Brillouin zone. 
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Figure 3: Low-lying 3D band structures for: (a) zigzag AA bilayer near the M point, (b) 
armchair AA bilayer near the Gamma point, (c) zigzag AAA trilayer near the K point, and 
(d) zigzag ABA trilayer near the K point. 
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Figure 4: The charge density p of (a) pristine AA stacking, and (b) zizag AA bilayer. The 
charge density difference A p of zigzag AA bilayer in (c) yz-plane, and (d) xz-plane. 


27 











DOS(states/eV.unit cell) DOS(states/eV.unit cell) DOS(states/eV.unit cell) DOS(states/eV.unit cell) 





Figure 5: DOS of the oxygen-absorbed few-layer graphenes: (a) zigzag monolayer, (b) 
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zigzag AA bilayer, (c) armchair AA bilayer, (d) chiral AA bilayer, (e) zigzag AB bilayer, 






































































